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Abstract 

In this paper we present the techniques for computing cosmological bounces in 
polynomial f{R) theories, whose order is more than two, for spatially flat FLRW 
spacetime. In these cases the conformally connected Einstein frame shows up multi¬ 
ple scalar potentials predicting various possibilities of cosmological evolution in the 
Jordan frame where the f{R) theory lives. We present a reasonable way in which 
one can associate the various possible potentials in the Einstein frame, for cubic 
f{R) gravity, to the cosmological development in the Jordan frame. The issue con¬ 
cerning the energy conditions in f{R) theories is presented. We also point out the 
very important relationships between the conformal transformations connecting the 
Jordan frame and the Einstein frame and the various instabilities of f{R) theory. 
All the calculations are done for cubic f{R) gravity but we hope the results are 
sufficiently general for higher order polynomial gravity. 


1 Introduction 

Due to advancement in technology we can presently probe most of the information stored 
in the cosmic microwave background radiation (CMBR) held. The average temperature of 
CMBR and the huctuations on it gives us a hint of the pre-CMBR phase of the universe. 
The data stored in CMBR hints at the possibility that the very early phase of the universe 
may have gone through an inhationary phase. The theory of cosmic inhation BE] tries to 
address various problems in cosmology, as the horizon problem, the entropy problem, the 
hatness problem and others. In spite of the fact that inhationary theories can solve most 
of the problems of the very early universe they have an inherent problem, termed as the 
trans-Planckian problem misiEiii]. To address the difficulties of the trans-Planckian 
problem and the major problem of the big-bang singularity there are proposals for an 
alternative way to deal with the very early universe. The alternative paradigm about the 
cosmological dynamics of the early universe deals with the idea of a cosmological bounce. 
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In the bouncing universe scenario there exists a contracting phase of the universe, prior 
to the time of the big-bang singularity. The contracting universe never actually reaches 
the singular point but bounces to an expanding phase. Technically the scale factor of 
the FLRW spacetime never vanishes but it attains its minimum value during the bounce. 
Some authors have proposed cosmological bounces in the general relativistic framework 
and others have used theories using modified gravity/particle physics to model such 
cosmological bounces [ini[IIlll2lll3l[Ill[l5lll6l[I7l[l8llil|20l|2T]. 

Cosmological bounces can be thought of as an interesting alternative to inflation as it 
resolves the horizon problem by making most part of the observable universe to be causally 
connected during the contracting phase of the universe. The perturbation length scales 
corresponding to these observations became superhorizon during the contracting phase. 
These perturbations later entered the causally connected universe in the expanding phase. 
One may assume that far away from the bouncing regime the dynamics of the cosmos is 
guided by general relativity (GR) and only when the Ricci scalar attains a considerable 
value (near bounce) the theory of gravity attains quantum corrections [22] and an effective 
f{R) theory of gravity emerges One can use f{R) theory of gravity near the bouncing 
regime to understand the dynamics near the bounce. It must be noted that in f{R) 
gravity theories one can have non-zero, and presumably, high values of the Ricci scalar, 
R, in a primarily radiation dominated (early) phase of the universe unlike GR, where the 
Ricci scalar is zero in radiation domination. The reason being that the dynamic equation 
for evolution of the metric changes in f{R) theories (the standard Einstein equation of GR 
is modified/altered). Due to this fact one can indeed think of high Ricci scalar values near 
a radiation dominated bouncing regime. Some preliminary work on bouncing cosmologies 
using f{R) theories can be found in Ref. |2Sl I2ZI [2S|- In a previous paper Ref. [22] the 
present authors proposed an interesting cosmological bounce in quadratic f{R) theory in 
presence of hydrodynamic matter. In the previous work it was shown that for spatially 
flat FLRW spacetimes one cannot get cosmological bounces in both the Jordan frame, 
where the f{R) theory lives, and the conformally connected Einstein frame. In spite of 
this fact one can use the Einstein frame to solve the bouncing problem. In the present 
work we extend our earlier work by touching upon the point related to energy conditions. 
We show that one can always use the Einstein frame to calculate the bouncing dynamics 
in general polynomial gravity without breaking the weak energy condition and the null 
energy conditions. It is also pointed out that some energy conditions as the strong energy 
condition and the dominant energy condition may be violated near the bouncing point. 
Actually the cosmological bounce only happens in the Jordan frame. In this paper when 
we refer to the bouncing point in the Einstein frame we mean the cosmological coordinates 
in the Einstein frame which corresponds to the cosmological variables in the Jordan frame 
at the time of bounce. The energy conditions are perfectly dehned in the Einstein frame 
where the mathematical description of the cosmological process follows the known route 
of GR. In the Jordan frame the energy conditions are violated but this violation may not 
be taken very seriously as the energy conditions are not uniquely formulated for f{R) 
theories as the energy momentum tensor gets contribution from curvature terms. In this 
paper we explicitly deal with spatially flat FLRW spacetimes where the dynamics of the 
scale factor is guided by cubic f{R) theory in the Jordan frame. The interested reader 
can see Ref. [3D] where the peculiarities of the two conformal frames are discussed. 

^Most of the common techniques of f{R) theories which we will use in this article are appropriately 
summed up in the review articles in, Ref. [23], Ref. [24] and Ref. [25] . 
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In this paper we address the important question relating to the conformal correspon¬ 
dence of the Jordan frame and Einstein frame where the form of f{R) is a polynomial 
function, whose order is more than two. This question is particularly important in f{R) 
theories as because the correspondence between the two conformal frames for such higher 
order f{R) theories becomes many-to-one. While there is only one f{R) theory in the 
Jordan frame, for a polynomial gravity theory, there can be various different possibilities 
of cosmological evolution in the Einstein frame. It appears that there are multiple Ein¬ 
stein frame descriptions of a single f{R) theory. This fact emerges from the conformal 
transformations connecting the two frames. For quadratic gravity this correspondence 
was unique and there was no confusion in solving the bouncing problem in the Einstein 
frame and then converting the results into Jordan frame language. But as soon as one 
starts with cubic gravity the problems of the conformal frames becomes apparent. In 
this work we address the problem of the emergence of multiple scalar potentials in the 
Einstein frame for a single polynomial gravity theory and specify how one can identify 
the origin of the multiple potentials in the Einstein frame. The discussion regarding mul¬ 
tiple Einstein frame description, of a unique f{R) theory, is general in nature and has no 
particular connection to bouncing cosmologies. We hope our results will give a new way 
of interpreting polynomial f{R) theories. 

In this paper we point out that there is a close relationship between the instabilities 
of f{R) theories and the conformal transformations connecting the Jordan frame and the 
Einstein frame. It is well known that if f'{R) < 0 the gravitational theory becomes diffi¬ 
cult to handle as the effective gravitational constant in f{R) theory diverges or becomes 
negative. On the other hand if f'{R) < 0 the conformal transformations connecting the 
two conformal frames becomes ill defined. Going one step further, it was shown by Dolgov 
and Kawasaki in Ref. ina that for low curvature f{R) gravity the condition /"(i?) < 0 
sets up a new form of instability where the effective gravitational interactions increases 
without an upper bound. Dolgov and Kawasaki’s result deals with low curvature cases 
and consequently it was not clear whether their results hold for f{R) theories applied in 
the very early universe where the Ricci scalar may not have a small value. In general 
in Ref. [30] the authors also proclaim that f”{R) > 0 is a stability condition for their 
cosmological model. In this article we show that in the case of the very early universe, 
where the gravitational dynamics is dictated by a cubic theory of gravity, the epoch when 
R attains the value Rc, where Rc is the solution of f"{R) = 0, is not connected with 
other cosmological epochs. In other words a bouncing universe will take an inhnite time 
to attain R = Rc after the bounce. Consequently the very early universe can exist in two 
complimentary branches, in one branch f"{R) > 0 and in the other branch f"{R) < 0. In 
the branch where f"{R) < 0 one can have perfectly well behaved cosmological bounces but 
this branch also shows the embryo of gravitational instability, as in this branch the Ricci 
scalar R can become arbitrarily small and consequently Dolgov-Kawasaki like instabilities 
can affect the theory in the low curvature limit0 On the other hand as we are applying 
f{R) theory to study cosmological bounce, which can only happen in the very early uni¬ 
verse and at high values of the Ricci scalar, the necessity of such a modified theory of 
gravity may decrease in the low curvature limit. Consequently near the neighborhood of 
R = 0 the theory of gravity may cross-over to conventional GR. If one uses cubic gravity 

^This branch will be governed by the scalar potential V 2 {(j)) in the Einstein frame as shown later in 
this article and this branch includes the point i? = 0 where R is the Jordan frame Ricci scalar. 
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in the other branch]^ then this theory does not have any apparent instabilities but the 
lowest value of Ricci scalar attainable in this branch is not zero and consequently it may 
be difficult to transform this effective f{R) theory into a conventional theory of gravity 
as GR. The two complimentary branches of cosmological existence meets at i? = Re- It is 
interesting to note that the Einstein frame description of f{R) gravity also breaks down 
precisely at the point R = Re- In this article we provide some general analysis on the 
relationship of the conformal transformations and the instabilities which plague any f{R) 
theory of gravity. As the instabilities are inherently related with the signs and zeros of 
the hrst and second derivative of f{R), with respect to R, we will analyze the properties 
of these derivatives and show how these properties affect the conformal transformations. 

The material in the present paper are organized as follows. In the next section we 
briefly recapitulate the conformal relationship between the Jordan frame and the Einstein 
frame for general f{R) theories. In this section we present the basics of the bouncing 
cosmological problem in f{R) theory and opine on the very important issues related 
to the signs and zeros of the hrst and second order derivatives (with respect to R) of 
f{R) theories and their relationship with the conformal transformations. The second 
section ends with a detailed discussion on the energy conditions in f{R) theories keeping 
the discussion focussed around cosmological bounces. The next section addresses the 
problem of multiple scalar potentials which arise in the Einstein frame when f{R) is a 
polynomial of order more than two. In this section we choose to work with cubic gravity 
where the problem of multiple scalar potentials hrst shows its signature. In this same 
section the analysis of the cosmological dynamics of f{R) theories around R = Rc is 
presented. In section the ways to understand the multiple potentials are discussed and 
explicit numerical solutions giving cosmological bounces in cubic gravity are presented. 
This section ends with a detailed discussion of the bouncing phenomena in cubic gravity. 
The last section concludes the paper. 

2 Metric f{R) gravity and its Einstein frame descrip¬ 
tion 

In this section we briehy describe the relationship of the two conformally related frames 
important for us. One is traditionally referred as the Jordan frame where the metric 
f{R) theory is dehned and the other is the Einstein frame where one can express the 
gravitational evolution of the cosmological system dehned in the Jordan frame. In this 
article we will be particularly interested in hat FLRW spacetimes and consequently the 
discussions will be based on equations of cosmological evolution where the curvature 
constant k has been set to zero. The elaborate relation between these two frames were 
described in an earlier paper Ref. [22] • Here we briehy express the most important results 
which connect these two frames. 

Using the FLRW spacetime in terms of the scale-factor a(t), 

ds^ = —dt^ -I- a?‘{t) \dr‘^ + r‘^{d6^ + sin^ 6 d(j)^)] , (1) 

^This branch is governed by the potential Vi{(j)) in the Einstein frame. 
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( 2 ) 

(3) 


one can write the field eqnations in metric f{R) gravity as |24] : 


3H^ + 2H 




— K 

F{R) 


PeS 5 


where H is the conventional Hubble parameter defined as H = a/a and k stands for the 
constant specifying the curvature of the 3-dimensional spatial hypersurface. The constant 
K = SttG where G is the universal gravitational constant. In the above equations 


F(R) 


df{R) 

dR 


(4) 


The dot specifies a derivative with respect to cosmological time t. The effective energy 
density, peff, and pressure, Pes, are defined as: 


PeS — P T Pcurv 5 PeS — P Pcurv ) 


(5) 


where pcurv and Pcurv are given by 


Pcurv 


Pc, 


RF-f_ 3HRF'{R) 

2k k ’ 

iPF" + 2HRF' + RF' RF - f 
K 2k ’ 


( 6 ) 

(7) 


which are curvature induced energy-density and pressure. In the above equations the 
primes designate a derivative with respect to the Ricci scalar R. The curvature induced 
thermodynamic variables exists in absence of any hydrodynamic matter, in contrast to 
the conventional p and P in 


= (p + P)u^Uu + Pg^u , (8) 

which has the information of hydrodynamic matter. In this article we assume the fluid to 
be barotropic so that its equation of state is 

P = up, (9) 


where a; is a constant and its value is zero for dust and one-third for radiation. It must 
be noted that in Eq. ([^ is the four-velocity of a fluid element and = — 1. Till 
now the description of the gravitational field equation and energy momentum tensor was 
specified in the Jordan frame. 

To understand the dynamics of f{R) gravity one can recast the problem in the Einstein 
frame by applying a conformal transformation on the Jordan frame metric as 


g^iu R{R)gfj,w , 


and simultaneously defining a new scalar field 0 as 


0 = 



InP(P). 


( 10 ) 


( 11 ) 
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This scalar field plays an important role in the Einstein frame. The conformally trans¬ 
formed line element in the Einstein frame is 


= —dP -I- a^(f) [dP -I- P{d9‘^ -I- sin^ 9 , (12) 

where the time coordinate, f, and the scale factor, d, in the Einstein frame are related 
to their corresponding Jordan frame terms via the relations di = ^JF{R) dt and d{t) = 
F{R) a{t). Using these transformations one can formnlate the gravitational dynamics 
of f{R) gravity in the Einstein frame in presence of matter and the scalar field 0 acting 
as sonrces. The energy-momentum tensor in the Einstein frame, which is related to 
in the Jordan frame, turns out to be 

T^u = ip + P)u^Uu + Pg^Lu , (13) 

where p = p/F‘^{R), P = P/F‘^{R) and = ■\/F{R)u^. In the Einstein frame = 

— 1. Except the energy-momentum tensor for the scalar field also acts as source of 
curvature in the Einstein frame and it is given as 


- jo£(0), 


(14) 


where the scalar field Lagrangian is 

^4>) = ^9<,<t>804‘3“^ + V{4i). 


(15) 


The scalar field potential in the Einstein frame turns out to be 


V(4>) 


RF-f 

2kF^ 


(16) 


where one has to express R = R{4>), from Eq. (0 by inverting it, and then express U(0) 
as an explicit function of 0. From the form of and the signature of the metric used, 
= —p and S\ = P yielding: 


P(j> 


1 / d0\ 

2\d!) 


+ V{4>), Pt=\ 



(17) 


where the scalar field 0 is assumed to be a function of time only. The total energy- 
momentum tensor responsible for gravitational effects in the Einstein frame is 
which is a mixed tensor with only diagonal components. 

The time coordinate, t, and the Hubble parameter, FI, in the the Jordan frame are 
related to the time coordinate, i, and Hubble parameter H{= 1^), in the Einstein frame 
via the relations: 


h = (18) 

As we will be interested mainly in bouncing cosmologies R will be set to zero. The 
instant tg = 0 is the bouncing time in the Jordan frame. The Einstein frame description 
of cosmology can be tackled like FLRW spacetime in presence of a fluid and a scalar 
field. The presence of the Scalar field potential V (0) gives one a pictorial understanding 
of the physical system which is lacking in the Jordan frame. Seeing the nature of the 
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potential and the initial conditions of the problem one gets a hint about the possible time 
development of the system. The time evolution of the scalar held in the Einstein frame 
is dictated by the equation 


, nfrd(j) , dV 

—X- + 6U -p H—— 

dt^ dt d(j) 


(1 -3ui)p, 


(19) 


where the equation of state of the huid in the Jordan frame is P = up. It is interesting to 
note that the equation of state of the huid remains the same in the Einstein frame. The 
evolution of the energy density in the Einstein frame is given by 




( 20 ) 


The above two equations dictate the time evolution of p and 0 in the Einstein frame. To 
generate proper bouncing solution from the above two equations one requires the values 
of only two quantities at the bouncing time, they are 0(to), other parameters 

are determined from these two at the bouncing time[29]. In general for a hat FLRW 
spacetime these two values of the respective quantities are enough to solve the whole 
system in the Einstein frame. The expression of the Hubble parameter and its rate of 
change in the Einstein frame are given as 


dH 

di 


^(pp + p). 



( 21 ) 

( 22 ) 


In the next section we will use the above results to formulate a version of energy conditions 
which can be applied in metric f{R) gravity. 

It is to be noted that the prescription for transforming from the Jordan frame to the 
Einstein frame is limited and fails in a particular case. In this paper we will show an 
example of this failure. More over the two conformal frames may not be equivalent for 
all cosmological scenarios, specially for bouncing cosmologies in hat FLRW cases in the 
Jordan frame there corresponds no bounce in the Einstein frame |29j . 


2.1 Cosmological bounce in Jordan frame for polynomial f{R) 

As in this paper we will be mainly talking about cosmological bounces so it is pertinent 
to give the conditions for a bounce. The conditions for a cosmological bounce are: 

P(fo) = 0, P(fo)>0, (23) 


where H is the Hubble parameter in the Jordan frame where the gravitational ehects are 
dictated by some /(P) theory of gravity. We will assume t = to to be the bouncing time. 
In this paper most of the time we will set the origin of cosmic time at the bouncing point 
and consequently to = 0. From Eq. (18) it can be seen that t = 0 in the Einstein frame 
at the time of bounce, t = 0, in the Jordan frame. 

The condition that the Hubble parameter vanishes at the bouncing time transforms 
to the condition: 


Po 


Po/o - fo 

2k, 


= 0 . 


(24) 


7 






for a flat FLRW spacetime. In the above equation we have used subscripts zero to de¬ 
note the values of the quantities at the bouncing time when t = 0. If one assumes the 
hydrodynamic matter in the cosmological background satishes the condition^ 

P>0, p + P>0, (25) 

then the other bouncing condition becomes 

^o/o' + ^o/o <0, (26) 

for the flat FLRW solution. In Ref. [22] it was shown that if one takes f{R) = R + aR? 

then the bouncing conditions invariably constrains a < 0 and a successful cosmological 

bounce only happens in the presence of hydrodynamic matter. 

In this paper we will focus on cubic gravity where 

f{R) = R + I3R^ + -iR\ (27) 


where f3 and 7 are real numbers. Using this form of f{R) one can easily show that there 
can be two kinds of bounces, one without any presence of matter and the other is a 
bounce in presence of hydrodynamic matter. In particular for the flat FLRW spacetime 
in absence of any matter the bouncing condition becomes 


(flo/i-/o) = 0, 


(28) 


predicting that 

= (29) 

27 

As Rq is positive dehnite at the bouncing point for flat FLRW spacetimes, one concludes 
that 7 and f3 must have different signs for cosmological bounces taking place in the 
absence of any hydrodynamic matter. In the present case if one demands that f'{R) > 0, 
the coefficients (3 and 7 must satisfy the following inequality 


0 < < 37. 


(30) 


From Eq. (29) it was seen that [3 and 7 should have opposite sign whereas from the above 
equation it can be uniquely said that for a bouncing solution 7 > 0 and /3 < 0 . 


2.2 A discussion on the signs and zeros of the first and sec¬ 
ond order derivatives of f{R), with respect to R, and their 
relationship with the conformal transformations. 

In this article we will be analyzing bouncing phenomena with polynomial f{R)- Polyno¬ 
mial f{R) theory themselves have some interesting features. To discuss those features we 

^It is to be noted that the conditions on the energy density and pressure assumed here does not stem 
from any energy condition and one cannot strictly call it the weak energy condition as the authors did in 
[29]. Although the conditions can be derived by assuming > 0 in the Jordan frame where is 

any time-like 4-vector, but then also we cannot strictly call this conditions the weak energy condition as 
Rjy is not the only source of curvature, as implicitly assumed in Einstein gravity, as there are effective 
Ricci scalar dependent energy density and pressure. The only way to justify the conditions is by assuming 
no exotic hydrodynamic matter during the bounce phase. Discussions on the energy conditions in f{R) 
gravity will be presented later in this paper. 




first state the basic issues about stability of a f{R) gravity theory. From Eqs. ([^ and 
(|^ it is seen that k,/F(R) acts like an effective gravitational constant in f{R) theory and 
consequently for F{R) > 0 the theory is well dehned. On the other hand if F{R) < 0 
the theory of gravitation becomes either ill dehned or produces a negative effective grav¬ 
itational constant which makes the theory unstable. On the other hand if F{R) > 0 
but F'{R) < 0 there can be other instabilities in the theory in the low curvature limit. 
This second kind of instability is generally called a Dolgov-Kawasaki instability |3l] . The 
authors who hrst analyzed this kind of instability used a weak gravitational held and 
applied f{R) theory where as in our case we will never use a weak gravitational held 
approximation. Consequently, we do not expect Dolgov-Kawasaki like instability to ahect 
the stability of the cosmological bounce phenomena. Unlike the Dolgov-Kawasaki case 
for low curvature gravity where one obtains unstable solutions for F'{R) < 0, in the early 
universe one can have both F'{R) > 0 as well as F'{R) < 0, both the branches giving 
rise to perfectly well behaved bounces. It will be shown later that in higher derivative 
cosmology guided by f{R) theory the universe avoids to attain the value of R which is 
a solution of F'{R) = 0. As a consequence the two branches are complimentary to each 
other. If the universe is in the F'{R) > 0 branch it will remain so for inhnite time and if 
it is on the other branch it will remain there for inhnite time. 

In f{R) theory the signs and zeros, of the hrst and second derivatives of f{R), plays 
an important part in the cosmological dynamics in the Jordan frame and the conformal 
transformations connecting these two frames. The ehect of the sign of the hrst derivative 
of f{R) has been discussed in the last paragraph, henceforth we discuss the ehects of 
the sign of the second derivative of f{R) on the conformal transformation. More over it 
will be seen that the roots of F'{R) = 0 plays an important role as far as the conformal 
transformations are concerned. Here we present three general results about polynomial 
f{R) theory and the conformal transformations relating the Jordan frame and the Einstein 
frame. Henceforth in this subsection we will assume f{R) to be a polynomial of order n so 
that F{R) is a polynomial of order n — 1 and F'{R) is a polynomial of order n — 2. The hrst 
point about an ehective Einstein frame description of an f{R) theory is as follows. Any 
Einstein frame description of gravity arising from an odd order polynomial f{R) where 
the order of the polynomial n > 1 must have multiple Einstein frame potentials. To prove 
this statement we hrst notice that if one wants to write an unique potential U(0), as 
given by Eq. (16), for the scalar held in the Einstein frame one has to invert the relation 
which gives 0 as a function of R in Eq. (11) and express R = R{4>). On the other hand 
it is seen from Eq. (11) that (j){R) = \J ^ \nF(R) is not uniquely invertible if F(R) is a 
polynomial of degree greater than or equal to two. It will be shown later that if F{R) is 
an odd order polynomial (hence f{R) is an even order polynomial) one can circumvent 
this problem related to invertibility. Only when F{R) is an even order polynomial (where 
f{R) is an odd order polynomial) the problem related to invertibility becomes really an 
acute one. If F{R) is a polynomial of degree greater than two then the inverse function 
R = R{4>) can be made single valued by partitioning it into various branches where in each 
branch the relation R = R{(j)) remains a single valued monotonic increasing/decreasing 
function of R. For all these separate branches there will be separate forms of U(0). If 
f{R) is an odd order polynomial, then F{R) is an even order polynomial in R and thus 
F{R) has at least one extremum for some hnite value of R. Consequently, the Einstein 
frame description for such an f{R) theory will have multiple scalar potentials U(0). If 
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Figure 1: Einstein frame scalar field as function of Jordan frame Ricci scalar. We have 
taken /3 = —10^^, 7 = in Planck units. 


F{R) has n extrema, then it has n + 1 different Einstein frame pictures separated by 
the points of extrema. If f{R) is an even order polynomial then F'{R) is also an even 
order polynomial which may not have any real zeros if the parameters of the theory are 
chosen properljj^ In such a case the question of multiple potentials may be avoided. In 
Fig.|I] we present the graph of the scalar field 0 versus the Jordan frame Ricci scalar R 
for cubic f{R) gravity as given in Eq. (27). The graph clearly shows the multi valued 
nature of the inverse function. The minimum of (j) appears on the right hand side, near 
Rc = 5.0 X 10“^^ where 0 = 0c- Here Rc is the solution of the equation F'{R) = 0 
for cubic gravity. The inverse function is invertible in the right hand branch or the left 
hand branch of the minima where the branches bifurcate at the minima. In plotting the 
above graph we have used Planck units where (3 = —10^® GeV“^ = —where 
Mp = 10^® GeV is the approximate Planck energy. To represent the numerical values in 
a brief and compact way all the dimensional quantities are scaled by Mp in this paper. 
Henceforth when we give the values of dimensional quantities it will be assumed that we 
are using the Planck units. 

Our next observation regarding polynomial f{R) theory is as follows. One cannot 
simultaneously have both F{R) > 0 and F'{R) > 0 for all possible values of R in a 
polynomial f{R) theory of gravity. To prove the above statement one must hrst note that 
if f{R) is an even order polynomial, then F{R) must be an odd order polynomial and 
consequently it must have at least one real zero. As a result F{R) > 0 does not hold for 
all R. On the other hand suppose f{R) is such an odd order polynomial so that the even 
order polynomial F{R) has no real roots and the condition F{R) > 0 holds for all R. But 
then F'{R) is an odd order polynomial and has at least one real root and consequently 
one cannot have F'{R) > 0 for all R. 


^In such cases there will be even number of roots of the equation F'{R) = 0 and one may in principle 
chose half of the roots to be complex numbers by tuning the parameters appearing in the algebraic 
equation. The other half of them will be complex conjugates of the earlier roots. 
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The third observation regarding the Einstein frame description of f{R) gravity follows. 
The Einstein frame description of f{R) gravity dynamics breaks down at the extrema of 
F{R). This observation about f{R) gravity is more general in nature and also holds 
for non-polynomial nature of f{R)- This observation was known previously in various 
forms and in Ref. [33] the authors gave a proof of the above statement. Here we present 
a different version of the proof which is particularly suited for f{R) theories governing 
cosmological dynamics. To prove it one must hrst note that the time evolution of the 
scalar held (j) in the Einstein frame is given by Eq. (19) which contains a term dV/dcf on 
the left hand side. One can calculate dV/dcf as 


dV{(j)) _ dV dR 
d(j) dR d(j) ’ 


where Eq. (16) can be used to calculate dV/dR and one can invert Eq. 0 to calculate 
dR/d(j). If it happens that F{R) has an extremum for some value of R then at that 
extremum point dcf/dR = 0 and a result dR/dcf is divergent. In general at these values 
of R where 0 has extremas, dV/dR are well behaved. Consequently dV/dcp becomes 
singular at the extremum point and the scalar held dynamics becomes ill dehned near the 
extremum point making the Einstein frame description of f{R) gravity invalid. 


2.3 A brief discussion on energy conditions in f{R) gravity and 
their conformal frame dependence 

The energy conditions on the energy-momentum tensor, for a perfect huid, are most 
well understood in Einsteins gravity where the energy-density, p, and pressure, P, of 
the hydrodynamic huid acts as source of curvature in space tim^ The curvature of 
space time back reacts on the source, via Einstein equation, and modihes them. As a 
consequence, in cosmology guided by GR, both the metric and the huid properties evolve 
in time. This evolution by itself can produce various forms of p and P which can produce 
peculiar situations where the energy density and pressure does not behave as they should 
for standard matter for which the energy conditions hold in GR. In such circumstances 
one may invoke exotic matter to justify the properties of p and P. The standard energy 
conditions in GR can be briehy described as [M] : 

1. The null energy condition (NEG) which is valid when p -|- P > 0. 

2. The weak energy condition (WEG) which is valid when p > 0 and p -|- P > 0. 

3. The strong energy condition (SEG) which is valid when p -|- P > 0 and p -|- 3P > 0. 

4. The dominant energy condition (DEG) which is valid when p > 0 and p > \p\. 

In the context of /(P) gravity there are various complexities in interpreting and applying 
energy conditions as the source of curvature of spacetime may itself contain the Ricci 

®In this section the variables p and P are used as generic symbols of energy-density and pressure 
in the Einstein frame. The same symbols were used earlier to specify energy-density and pressure in 
the Jordan frame. We are deliberately using these symbols as we do not want to increase the number 
of mathematical variables used in the paper and we think the reader can understand there role in the 
present context. 
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scalar R as seen in Eq. (|^ and Eq. (|^. In Ref. [35] the authors tried to implement the 
WEC and DEC in the Jordan frame using pefr and Pgff as defined in Eq. (|^, while the NEC 
and SEC were derived using the properties of Raychaudhuri’s equation. In Refs. |3Sl EZj 
the authors give a much vivid discussion on the properties of energy conditions in f{R) 
gravity. In these references the authors note that the energy conditions may not be 
simultaneously satished in both the Jordan and Einstein frames. 

On the other if one started with minimally coupled scalar held in the Einstein frame 
and then tried to conformally transform the theory in the Jordan frame then matter 
becomes non-minimally coupled to curvature in Jordan frame. In Ref. [3H| the authors 
tried to reformulate the NEC in the Jordan frame in such a way that it reduces to the 
normal NEC in the Einstein frame under appropriate conditions. 

In this subsection we will specihcally give an example where most of the energy con¬ 
ditions are violated in the Jordan frame, at the time of bounce, and none of the energy 
conditions are violated in the Einstein frame, at the corresponding time. 

The structures of Eq. (|^ and Eq. (|^ and the bouncing conditions in Eq. (23) imme¬ 
diately shows that for a successful bounce in the Jordan frame 

• one must have at the time of bounce (pefr)o = (Pcurv)o + Po = 0 where we have 


assumed p > 0 in Eq. (25). The subscript zero specifies the value of the quantities 
at the time of bounce. 

• The other condition at the bouncing point is (Pefr)o = Po + (Pcurv)o < 0. This 
implies (Pcurv)o < 0 as we have chosen Pq > 0. 

These conditions show that the curvature energy is at its lowest value during bounce, 
as the energy density of normal matter is maximum during the bounce. In general after 
bounce p gradually fades away but pcurv increases and drives the expansion of the universe. 
Consequently one can say that a cosmological bounce in presence of matter violates the 
energy conditions in the Jordan frame, if one tries to understand the energy conditions 
there in terms of the effective energy density and effective pressure. On the other hand for 
a bounce in the absence of hydrodynamic matter pcurv = 0 and Pcurv < 0 at the bouncing 
point. Consequently in this case also all the energy conditions are violated in the Jordan 
frame if one tries to understand the energy conditions there in terms of Pes and Pgfj. On 
the other hand in the Einstein frame the analysis of the energy conditions are different. 
We give the Einstein frame analysis of the energy conditions below. 

The correspondence between the /(P) theory variables in the Jordan frame in presence 
of a perfect fluid and the conformally related Einstein frame variables was presented 
before. As we have assumed normal hydrodynamic matter in the Jordan frame, whose 


properties were given in Eq. (25), it can easily be shown that in the Einstein frame one 
must also have 

p>0, p + P>0. (31) 

The above condition does not guarantee that the various energy conditions will be main¬ 
tained in the Einstein frame. The reason being that in the Einstein frame one has to 
apply the energy conditions on and consequently the quantities which go inside 

the energy conditions are pe = p + p</, and Re = P + P<f, where p^ and P^ are the energy 
density and pressure of the time dependent scalar held given as: 

2 1 / 7 ( \ 2 

- ^(0) • (32) 




p = I (^ 
2 V dt 
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From the above equations one can notice that + > 0 always in the Einstein frame but 

one cannot guarantee that > 0 in the Einstein frame and consequently the total energy 
density pe may turn out to be negative in the Einstein frame. This can happen because 
Prf) consists of the potential E(0) as given in Eq. (16) which may attain arbitrary small 
negative values. If one imposes the condition that a meaningful description of Jordan 
frame cosmology can only be described by its Einstein frame description if V (0) > 0 then 
one can be sure that the WEC and NEC will be satished in the Einstein frame for most 
of the cosmological processes. For bouncing cosmologies, where the bounce happens in 
the Jordan frame, one cannot impose this condition in the Einstein frame, for polynomial 
f{R) gravity, as in these cases the Einstein frame potential generally is always negative for 
some values of 0. In spite of this fact one can work out the cosmological bounce problem 
in f{R) gravity using the Einstein frame, without breaking any energy conditions at the 
bouncing point. One can always work in the Einstein frame in such a way that the 
cosmological system respects all the energy conditions at the “bouncing point” (which 
corresponds to the actual bouncing point in the Jordan frame), but slightly away from 
the bouncing point some of the energy conditions, as the SEC or the DEC may be violated. 

To elaborate the above mentioned points we specihcally show the various kinds of 
possible bounces in the Jordan frame and their Einstein frame analogue. First let us 
concentrate bounce in the presence of matter. From Eq. (21) one can see that at the time 
of bounce, one must have 


v(4,) = -p. 


(33) 


As p > 0 we have p > 0 in the Einstein frame and consequently at the bouncing time, 
E(0) < 0 in the Einstein frame. If near the bounce point one can ensure p + > 

|I/( 0 )| then Pe > 0 , and as a result the energy conditions will be maintained in the 
Einstein frame and a successful bounce can happen in the Jordan frame where the energy 
conditions, interpreted in terms of the effective thermodynamic variables, are violated. 

For bounce in the absence of of matter, where p = p = 0, which is a pure curvature 
driven bounce we have from Eq.(33), the potential at the time of bounce as E(0) = 0[^ 
In this case at the time of bounce, in general 0(0) 7 ^ 0 in the Einstein frame. As bounce 
in absence of matter requires E( 0 ) = 0 at the time of bounce, one can note that the 
equation of state for such a scalar held in the Einstein frame will be 


UJ4, = 


which implies 


i0-^(O)-E(0) 
i0'2(O) + E(0) 


|^(^( 0 ) —)■ 1 , 


(34) 


which implies a kinetic energy driven free scalar held in the Einstein frame. This result 
is in general true for all matter less bounces, which can be tackled in the Einstein frame, 
for spatially hat FLRW spacetimes. More over as in this case E(0) = 0 at the time of 
bounce f = 0 , 

. 


^Although a cosmological bounce, in the absence of any hydrodynamic matter, in the Jordan frame 
requires V (0) = 0 but the vanishing of the Einstein frame potential does not necessarily imply a cosmo¬ 
logical bounce, in the absence of hydrodynamic matter, in the Jordan frame in general. 
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If one puts the condition that 0^0) = 0 then it should have been a symmetric bounce in the 
Jordan frame but the peculiarity of cosmological bounces, in absence of hydrodynamic 
matter, makes these symmetrical bounces in the Jordan frame very rare. The reason 
being that if 0^0) = 0 then = |0'^(O) + ld(0) < 0 just near the bounce point as 
there 0'(O) —)■ 0 and V{(j)) < 0 breaking the energy conditions in the Einstein frame, 
which makes the dynamic development of the system impossible. Consequently any pure 
curvature driven bounce, where energy conditions are maintained in the Einstein frame, 
mostly are asymmetric in nature in the the Jordan frame. 

From the expressions of p^j, and it is clear that if V (0) < 0 and consequently 

DEC may be broken near the bounce point, in the Einstein frame for pure curvature 
driven bounce. On the other hand may be negative if E(0) > 0 and consequently 
in such regions SEC may be violated. It must be noted that the NEC and WEC are 
always maintained in the Einstein frame and none of the violations, of the other energy 
conditions, happen at the bouncing point, but may be violated near to it. 

The above examples show that one can have perfect cosmological bounces in the 
Jordan frame, for the spatially flat FLRW spacetimes, for which the analogous Einstein 
frame descriptions may not break any of the energy conditions at the bouncing point. 
The energy conditions in the Einstein frame may dictate the kind of bounce in the Jordan 
frame. 


3 The issue about multiple scalar potentials in the 
Einstein frame 

In this section we present discussions on the most difficult issue regarding the Einstein 
frame description of polynomial f{R) theories. The issue is related with the emergence of 
multiple scalar potentials E(0) in the Einstein frame for a single f{R) theory. To track 
the dynamics of f{R) gravity in the Einstein frame one has to solve the held equation for 
0 which requires the knowledge of V (0). From Eq. it is seen that if one can construct 
a potential E(0) in the Einstein frame using Eq. (16) then one can uniquely determine 


the behavior of 0 in the Einstein frame and once this is done one can transform the results 
back to the Jordan frame. This program becomes difficult if one is unable to hnd out an 
unique scalar potential in the Einstein frame corresponding to a unique f{R) theory in 
the Jordan frame. In this section we will hrst discuss how multiple potentials arise, in 
the Einstein frame, in the case of cubic f{R) theory of gravity. After pointing out the 
difficulty posed by this multiple Einstein frame description of polynomial f{R) theories 
we will try to address how these difficulties can be addressed. 

To illustrate the difficulties inherently related to polynomial f{R) theories we study 
the particular case of cubic gravity . In this case we have 

f{R) = R + ^R^ + -fR^, 


as given in Eq. (27). In this case one can easily check, from the relation 





that 


= 0 , 

So' 3o' V / 
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Figure 2: The nature of the approximate unique Einstein frame potential for f3 = —10^^ 
and 7 = in Planck units.This potential only reflects the true dynamics of cubic f{R) 
theory when i? — 0 . 


gives R in terms of 0. This relation is not linear and yields two roots as 


- - I V 




(35) 


which shows where Ri > R 2 where the plus sign goes with the hrst branch. For a 
cosmological bounce in the Jordan frame it was discussed before that 7 > 0 and (3 < 0. 
The above results show that R = R{4>) is not a single valued function. One can tackle this 
problem by using two branches for R where we have for branch one Ri = Ri{4>) and for 
the second branch R 2 = i? 2 ( 0 )- Each of these branches are single valued functions of 0. 
It is interesting to note that on the hrst branch f"{Ri) = F'{Ri) > 0 and for the second 
branch /"(i? 2 ) = F'{R 2 ) < 0. These two branches meet at a particular 0 called 0c whose 
value will be deduced shortly. These two branches gives two different possible values of 


V (0) in the Einstein frame if one uses Eq. (16) to express the form of the potential in the 
Einstein frame. 

In the case of cubic gravity 


RF - f = (3R^ + 2jR^ 


There is an interesting approximation using which one may avoid the problems arising 
out of the multiple possible scalar potentials in the Einstein frame. If one just neglects 
the cubic term in R while inverting the relation between 0 and R then this gives rise to 
the approximate but unique result 

20 V J 

Now one can write an unique scalar potential for cubic gravity in the Einstein frame as 


1 /( 0 ) ~ 


2k [40 V ) 


403 V 


(36) 
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Figure 3: The two branches of the exact Einstein frame scalar held potential for /3 = 
— 10^^, 7 = in Planck units. The solid curve is V2(0) and the dashed curve is Vi(0). 
Here dVi^ 2 /d(j) —> oo a.t (f) = 4>c- The potentials do not exist for 0 < (pe¬ 


as obtained in Ref. [32] • From the form of the above potential one can verify that V (0) < 0 
for all 0 as /3 < 0 and 7 > 0. The above potential has only one zero for hnite values of 0 
and that zero is at 0 = 0. From now onwards we will call this Einstein frame potential 
as the approximate unique Einstein frame potential. This approximate unique Einstein 
frame potential is dehned for all values of 0 although it truly rehects reality only in the 
neighborhood of 0 = 0. The nature of the approximate unique Einstein frame potential 
for cubic gravity is shown in Fig. In this potential only the region near R = 0 is relevant 
and later we will show that this region also admits a cosmological bouncing solution in 
the corresponding Jordan frame. 

On the other hand if no approximations are made then one has two scalar potentials 
1 / 1 ( 0 ) and 1 / 2 ( 0 ) in the Einstein frame corresponding to the cubic gravity Lagrangian in 
Jordan frame. The potentials can be written as 



where subscript 1 corresponds to the upper minus sign and 2 corresponds to the lower 
plus sign on the right hand side. These potentials correspond to the two branches of R 
in the R — cp plane. The above forms of the scalar held potentials in the Einstein frame 
immediately shows that as soon as 

, (38) 

the Einstein frame description of cubic gravity ceases to exist as the potentials turn out 
to be complex. From the expressions of Vi,2(0) it is clear that 10(0) < 0 for all values 
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of (j), for the bouncing solution, where (3 < 0 and 7 > 0. More over V2(0) has a zero at 
0 = 0. On the other hand Vl( 0) can have both positive and negative values for the above 
choice of parameters and it has a zero for 0 7^ 0. It can also be noted that both Vi(0) 
and 1/2(0) are zero as 0 —)■ cxd. 

As an exponential function cannot ever be negative so the term on the right hand side 


of the inequality, in Eq. (38), cannot be negative if the inequality holds. If one chooses the 
parameters 0 and 7 in such a way that the (1 — (0^/87)) is negative then the inequality 
never holds and the Einstein frame description of cubic gravity holds for all values of 0. 
If one wants to make the Einstein frame description of cubic gravity valid for all values 
of 0 then another issue related to the conformal transformations become problematic. 
The Einstein frame description of cubic gravity holds as long as F{R) > 0, otherwise 
one cannot dehne the scalar held itself. More over if F{R) < 0 then the original f{R) 
theory of gravity itself becomes unstable. If one requires F{R) > 0 then the inequality 


in Eq. (30) has to be satished and in that case the right hand side of the inequality in 


Eq. (38) become positive and consequently there arises a hnite value of 0 as 


2k \ 37 



(39) 


which sets the lower limit of the scalar held strength in the Einstein frame. The conformal 
transformations will never be able to generate any Einstein frame dynamics of cubic 
gravity if 0 < 0c. At this value of 0, which corresponds to 


Rr = -{ — 

37 


(40) 


one can verify that F'{R) = 0 and as a consequence the Einstein frame description of cubic 


gravity becomes ambiguous as discussed in subsection 2.2, The two potentials Vi(0) and 
1/2(0) are in general diherent from each other and Vi(0) = 1/2(0) only when 0 = 0c and 
0 —?• cx) i.e., at the two extremities of the range of 0. The approximate scalar potential 


1/(0), as given in Eq. (36), approximately equals 10(0) very close to 0 = 0. For this choice 
one can see that the Ricci scalar in the Jordan frame approaches zero as 0 —)■ 0 in the 
Einstein frame. Consequently the approximate potential, as proposed in Ref. |32] only 
holds near the i? —)■ 0 limit in the Jordan frame and for all other hnite values of the Ricci 
scalar one has to rely on the dynamics produced by one of the actual potentials 10,2(0) 
in the Einstein frame. The nature of the two potentials 10(0) and 10(0) for cubic f{R) 
gravity are shown in Fig. The plots show that the two potentials meet at 0c which is 
negative for the particular choice of parameters 0 and 7 in cubic gravity. The potentials 
also merge at 0 —>■ 00 which can be inferred from the plots as the two potentials converge 
as 0 increases in the positive side. The hgure shows that 10(0) is always negative except 
at 0 = 0 and 10(0) has a zero for a non-zero 0. 

The above arguments show that the Einstein frame description of cubic gravity (de- 
hned in Jordan frame) is limited. On the other hand one may partition the Einstein 
frame description of cubic gravity in two regions dictated by 10(0) and 10(0) and ap¬ 
ply the conformal transformations safely but as soon as one reaches 0c, where these two 
potentials merge with each other, the Einstein frame description fails. The next section 
shows that the Jordan frame analysis of the cosmological system near R = R^, where the 
Einstein frame description fails, reveals an interesting feature of cubic f{R) theory. Using 
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the approximate potential, as given in Eq. (36), one never feels this constrained nature 
of cosmological dynamics in the Einstein frame, as in this approximation the potential is 
valid for all values of 0, although it truly describes reality only near 0 = 0. 


3.1 The cosmological behavior of general polynomial f{R) the¬ 
ories very near to R = Rc in the Jordan frame 

We have seen in the above discussions that the point R = Rc, where F\Rc) = 0, has some 
special properties as far as conformal transformations are concerned. In this subsection 
we will elaborate a bit more about the specihc nature of cosmological dynamics, in the 
Jordan frame, when the system is near R = Rc- This analysis is necessary as near R = Rc 
the conformal transformations fail and one cannot generate any Einstein frame analog of 
the Jordan frame dynamics and consequently one has to rely only on the Jordan frame 
dynamics of the system. For this analysis we will employ linear perturbation theory near 
R = Rc where the Hubble parameter is He- As Eq. ([^ and Eq. ([^ are written purely 
in terms of R, H and p, for the spatially flat FLRW case, it is enough to specify these 
parameters near {Rc, He, Pc) point, in the R — H — p space, as 

R = + SR ) 

H = Hc + 6H, 

p = Pc + Sp 


where SR , SH and Sp are small fluctuations of R, H and p over Rc, He and pc for any 
arbitrary f{R) theory. To linear order in fluctuations one can write 

f{R) = f{Rc) + F{Rc)SR, 

near the {Rc,Hc) point where F{Rc) = f'{R)\R=R^- As we have assumed F'{Rc) = 0 we 
have, 

F{R) = F{Rc), 

F\R) = F"{Rc)SR, 

F"{R) = F"{Rc) + F'"{Rc)SR, 


near the {Rc,Hc) point. In this analysis we will assume that F{R) > 0 as this ensures 
stability of the f{R) theory. Using the above results in Eq. (|^ and Eq. ([^ and calculating 
order by order one obtains 


and 


3H^c 

HeSH 

H 


K 

F{Rc) 

K 

QF{Rc) 

SH = - 


Pe+j^{RcF{Rc)-f{Rc)} , 
Sp, 


+ u) 

2F{Rc) 


Sp, 


(41) 

(42) 

(43) 


Pc(l + <^) — 0 . 


(44) 
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Figure 4: Bounce with radiation background in branch V2(0), shown by the dashed curve, 
and in branch 4h(0), shown by the solid curve. In these plots {3 = —10^^, 7 = in 
Planck units. 


The above set of equations reveals more information if one rewrites Eq. (41 )in the following 
form: 

K 


3H3 = 


F{Rc) 


[p, + F'‘(K)V (,!,,)] 


(45) 


where in the above equation we have used the expression of 14(0), as given in Eq. (16), an 
Einstein frame object deliberately used for a specihc reason. For cosmological bouncing 
scenarios in the absence of any hydrodynamic matter, in cubic gravity, one has V (0c) < 0. 
As a consequence the above equation predicts that the universe will never be able to reach 
the R = Rc point after bounce as at that particular cosmological epoch < 0. The 
discussion below will show that not only for bouncing solutions in vacua, or any other 
solutions with hydrodynamic matter, where Eq. (45) gives > 0, the evolving universe 
will never be able to reach the phase when the Ricci scalar attains the value Rc- 

To understand the general nature of cosmological evolution near R = Rc one has to 
concentrate on the physical content of Eq. (44). This equation states that very near to 
the {Rc-, He-, pc) point either the universe is devoid of hydrodynamic matter or the matter 
present must be acting like a fluid with an equation of state a; = — 1. These two cases 
can be separately analyzed. First we assume that the universe is hlled up with matter 
where pc 7^ 0 but a; = — 1. In this case one attains a local de Sitter universe around 
{Rc, He, Pc) point as now 5H = 0 from Eq. (43). Much away from the {Rc, He, Pc) point 
it is difficult to predict the nature of the Hubble parameter but one can always predict 
uniquely the value of u as it is a constant in our analysis. In this case it is observed that 
if the universe at all comes near the R = Rc point then that universe must have a local 
de Sitter phase. If one starts to track the cosmological evolution from some point away 
from the neighborhood of the {Rc, He, Pc) point one would see that the time required for 
such a development is inhnite. The reason for this observation can be understood by the 
equation governing time development which is given by 

dH 


t = 


H 
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As the universe nears the {Rc, He, Pc) point if —)■ 0 and consequently t —)■ cx). This implies 
that for cubic gravity the universe will take an inhnite time to reach the {Rc, He, Pc) point 
after bounce. 

If one particularly concentrates on a bouncing universe devoid of any matter, we have 
shown previously that the {Rc, He) point is never attained as at that point He turns out 
to be complex. In this particular case in between the bouncing point where H = 0 and 
if > 0 and the {Rc, He) point the universe has to pass through a phase where if = 0 as 
because H was increasing with time after the bounce and must have decreased to zero at 
some time after bounce so that H^ becomes smaller than zero near the {Rc,Hc) point. 
More over at the point where H = 0 one must have if 7^ 0 so that f? 7^ 0. This fact 
can be understood from the nature of the curve in Fig. [T| where a bounce in absence of 
hydrodynamic matter can happen for regions on the right hand side of the minima of 0 
(in the Vi(0) branch) and their f2 7^ 0. One can now apply the same logic around the 
point where ff becomes zero and show that the universe will take an inhnite time to reach 
that point. 

On the other hand if the universe had explicit hydrodynamic matter during bounce, 
whose equation of state ca 7^ — 1, then one can see from Eq. (44) that the energy density 
becomes exactly zero at the {Rc, He, Pc) point. This can only happen when the universe 
becomes inhnite in volume (or the physical volume which corresponds to a unit comoving 
volume diverges) such that the matter energy density becomes exactly zero. If we assume 
that at a hnite cosmological time after the bounce the physical size of the expanding 
universe to be hnite, then the {Rc, He, Pc) point can only be reached by the universe after 
an inhnite time. The density p, of hydrodynamic matter, will in general decrease as the 
universe expands after the bounce and tries to reach the {Rc, He, Pc) point, but it will 
never be exactly zero at any hnite time. Consequently in all of the cases one can observe 
that the bouncing cosmological system, guided by cubic gravity, will never reach the 
{Rc, He, Pc) point after a hnite cosmological time from the bounce. The {Rc, He, Pc) point 
in a certain sense will always remain inaccessible to the bouncing solutions in cubic gravity. 
The Einstein frame description also precisely fails at these points. The discussion of the 
nature of the cosmological system near the {Rc, He, Pc) point is interesting and we have 
probed cosmological evolutions near this point for bouncing solutions in cubic gravity. 
We expect our results to be more general and valid for other forms of polynomial gravity 
as well. 


4 Ways to choose between the three potentials in cu¬ 
bic gravity 

The natural question which arises in the case of polynomial f{R) gravity is related to the 
proper choice of the potential in the Einstein frame. The problem becomes non-trivial in 
the simplest of the cases in cubic gravity where one can have three different potentials in 
the Einstein frame which corresponds to a unique f{R) in the Jordan frame. The hrst 


amongst the three is the approximate unique Einstein potential as given in Eq. (36). This 
potential truly reflects cubic gravity for small values of the Ricci scalar R. The other 


two potentials are Vi(0) and V2(0) as given in Eq. (37). This potentials are in general 
valid for all values of R when either f"{R) > 0 or f"{R) < 0. In the following part of 
this section we show that there can be interesting physics involved in choosing the proper 
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Figure 5: Cosmological bounce in vacuum in branch Vi(0) where f3 = —10^^, 7 = in 
Planck units. 


potential in the Einstein frame. We will show that presence of matter during bounce can 
be an important factor which helps us to choose amongst the three possible potentials in 
the Einstein frame. 

In cubic gravity one can have a cosmological bounce in the Jordan frame in presence 
of matter. In this case one has to satisfy the condition given in Eq. (33). Following the 
discussion after Eq. (33) in subsection 2.3 it is seen that for a successful cosmological 
bounce in the Jordan frame one must have V{(f)) < 0 in the Einstein frame. As in this 
case a bounce always happens when V (0) < 0, we can infer that such a bouncing scenario 
may be described by both Fi(0), V2(0) or E(0) as given in Eq. (36) in the Einstein frame, 
as all of these potentials can have regions where V((/)) < 0. The behaviors of the scale 
factor for the flat FLRW spacetime near cosmological bounces in presence of radiation are 
shown in Fig. It is seen that both the potentials are capable of producing cosmological 
bounces in presence of matter. 

On the other hand if the universe is devoid of any hydrodynamic matter then from 
Eq.(33) one can see that at the time of bounce V((/>) = 0. In this case at the time of 
bounce, in general 0 7^ 0 in the Einstein frame. Consequently, pure curvature driven 
bounce must happen through the potential Vi(0) in the Einstein frame as only Vi(0) does 
have a zero for non-zero 0. In this case if dcj)/dt is non-zero when Fl( 0) = 0, the scalar 
held climbs to positive values of Fi(0) during the bouncing period. The plot showing the 
evolution of the scale factor a(t), for cubic gravity, near a cosmological bounce in vacuum 
is shown in Fig. 

From the above description of bouncing phenomenon in the Jordan frame we observe 
that although the approximate unique Einstein frame potential can describe a cosmolog¬ 
ical bounce in presence of matter but it has no information about pure curvature energy 
dominated bounce. Consequently one can infer that the scope of the approximate unique 
Einstein frame potential for describing cosmological bouncing phenomenon in the Jordan 
frame is very limited. It is only a true rehection of cubic gravity for very small values 
of the Ricci scalar R and more over it can only describe a matter induced bounce in the 
Jordan frame. The nature of a cosmological bounce in cubic f{R) theory where one em- 
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Figure 6 : Bounce with radiation background using the approximate unique Einstein frame 
potential with /3 = —10^^, 7 = in Planck units. 


ployes the approximate unique Einstein frame potential is shown in Fig. This bounce 
can only happen in presence of matter background and for our case we have chosen the 
background hydrodynamic fluid to be comprised of radiation. 

4.1 Detailed description of the cosmological bounces in cubic 
f{R) gravity 

Discussions regarding pre-big bang cosmology which discusses about cosmological bounce 
was noted long ago, in 1992, in the Ref. [39j. It is well known that the cosmological 
behavior can be quite different in the Jordan frame and Einstein frame. In particular, it 
was mentioned in [29] that a bounce in the Jordan frame is not usually associated with a 
corresponding bounce in the Einstein frame for flat FLRW cosmologies. In spite of these 
observations it is surprising that the Jordan frame bounce scenario can be understood 
from the Einstein frame scalar potential picture. 

The evolution of the universe near the bouncing point in the Jordan frame can be 
understood from the nature of Fig. [T] where the minimum of (j) occurs at R = Re- On the 
right of Rc in Fig. [^one sees that 0 is an increasing function of R. This can only happen for 
the branch where In F{R) is an increasing function of the Ricci scalar, and from the initial 
discussions in section]^ this happens in the branch Vi(0). Consequently the branch Vi(0) 
corresponds to that part of cosmological dynamics in the Jordan frame where R > Rc- 
In this branch R ranges from R = Rc to R = 00 . The other potential branch V2{4>) 
corresponds to the region R < Rc in Fig. where 0 is a decreasing function of the Ricci 
scalar. This branch does have the embryo of Dolgov-Kawasaki like instabilities when the 
Ricci scalar becomes small. On the other hand as because in this branch the Ricci scalar 
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can come close to zero, the original f{R) theory may cross-over to conventional GR in the 
neighborhood of i? = 0. From this discnssion one can trace the behavior of cosmological 
bonnce in the Jordan frame if one follows the behavior of the scalar field 0 in the Einstein 
frame. 

Let ns hrst discnss abont the particnlar potential branch Fi(0) in the Einstein frame, 
which can snpport cosmological bonnces in presence or absence of hydrodynamic matter. 
First we discnss abont cosmological bonnces in the absence of any hydrodynamic matter. 
Bonnce in the absence of matter in the Jordan frame occnrs when the scalar potential 
in the Einstein frame vanishes, as discnssed previonsly. In this branch, at the bonncing 
time t = 0, the valne of the scalar field 0 = 0o when Vi(0o) = 0. In this branch, there 
are several possible lines of development of the cosmological system. These possibilities 
depend on the “velocity” {d(j)/dt or 0') of the scalar field at the time of bonnce. The 
contraction phase can start from 0 < 0o, cross that point (Jordan frame bonnce), rise 
np to a maximnm along the potential Vi(0) and come down again. If the scalar field 
has enongh initial kinetic energy it can also cross over the maxima of Fi( 0 ) and then 
slowly roll towards 0 —)■ oo. If it chooses to come down again it crosses 0 = 0o again, 
bnt this does not correspond to a bonnce, which can be nnderstood by a carefnl look at 
the Einstein frame initial conditions. Then it goes down towards 0c. While the nniverse 
moves towards 0c in the Einstein frame, the cosmological system in the Jordan frame 
moves so that the Ricci scalar their becomes Re- In doing so the nniverse will reqnire an 
inhnite time in the Jordan frame. One can also program bonnce, in absence of matter, 
in snch a way that the Einstein frame scalar field starts from above 0 > 0o dnring the 
contracting phase, then the cosmological system crosses the bonnce point and goes down 
towards 0c. Cosmological bonnce in presence of hydrodynamic matter follows a similar 
conrse except for the fact that at the time of bonnce the scalar held assnmes some valne 
0 = 00 where Vi(0o) < 0. More over, one mnst remember that bonnce in the presence of 
matter can also ocenr in the V 2 ( 0 ) branch of the potential. 

Similar consideration holds for the potential branch V2(0) bnt only for matter indneed 
bonnce. In this case 0 assnmes some valne 0 = 0o at the time of bonnce, snch that 
h 2 ( 0 o) < 0. In this branch increasing or decreasing of 0 correspond to decreasing or 
increasing of the scalar cnrvatnre R. The particnlar case of bonnce in the V2(0) branch, 
where 0 starts from some 0 < 0o and evolves towards 0 = 0, is interesting. In this case 
the contracting phase may proceed to the bonnce point 0 = 0o and then crosses this 
point and goes towards 0 = 0. In this branch as i? —)■ 0 as 0 —?• 0, as apparent from Fig 
one can assnme that the nniverse is evolving to a radiation dominated phase where 
the dynamics is governed by general relativity. The region near the i? —)■ 0 can act as a 
cross-over region from f{R) theory to general relativistic dynamics. 

Before we end this discnssion we present the varions valnes of the scalar held and its 
time derivative nsed in the Einstein frame to nnmerically solve the eqnations dictating 
the bonncing phenomena. The plot in Fig. shows a cosmological bonnce in absence 
of hydrodynamic matter where 0o = —.11, 0o = 10“®]^ Note that ont of these two 
conditions only 0' can be chosen snitably. The valne of 0 at bonnce is already specihed 


by the relation in Eq. (29) which gives. 




( 46 ) 


Mhis implies 00 = -.05Mp = - .05 x GeV; Rq = = -lO^^ GeVh 
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In this case 0 starts from some point where Vi(0) < 0 (see Fig. [^, crosses cpo, rises up to 
a maximum along the potential and then comes down again. 

The plots in Fig. |^and Fig. shows cosmological bounces in presence of radiation in 
both the branches Fi(0), V2(0) and in the approximate unique Einstein frame potential 
for the same initial conditions 0o = —-13, 0 q = 10“®. Unlike the matterless case, here 
both the initial conditions can be chosen suitably, with the only condition that the scalar 
held potential at bounce has to be negative. In these cases 0 starts from some point 
(f) < —.13, crosses the point (j) = 0o, goes up to a maximum along the potential and 
then comes down again. The important thing to be noticed in the above mentioned 
bouncing phenomena, in presence of matter, is that the same conditions at the time of 
bounce produces three different kinds of time dependence of the scale factor a{t) in the 
Jordan frame. Consequently these bounces are distinctly different and shows the richness 
of polynomial f{R) gravity. 


5 Conclusion 

In this paper we have extended our analysis on cosmological bounces, for the spatially 
hat FLRW solutions, in higher derivative f{R) theories of gravity. In our earlier work [22] 
we proposed a scheme by which one can study the phenomena of cosmological bounces, 
using spatially hat FLRW solutions, in f{R) theories using only the equations of general 
relativity. The method relied on the validity of the conformal transformations which 
connected cosmological dynamics in the Jordan frame, where the dynamics is governed 
by f{R) theory, and the auxiliary Einstein frame which is used to do the calculations 
related to the problem. Once the problem is solved in the Einstein frame one can map 
the solutions to the Jordan frame and get the bounce dynamics. This program can be 
hawlessly applied for the case of quadratic gravity where f{R) = R + aR^ where a 
is some constant negative real number. In the case of quadratic gravity it turns out 
that there is only one unique Einstein frame description of the bouncing phenomena. 
It was immediately noticed that the methods used in the earlier publication were not 
adequate when one uses higher order polynomials as cubic f{R)- The Einstein frame 
picture of cubic gravity and other higher order polynomial gravity does not have a one- 
to-one correspondence with the Jordan frame description of the cosmological evolution and 
consequently one does not know a priory which picture in the Einstein frame corresponds 
to the actual picture in the Jordan frame. In the present work we have proposed a 
consistent method using which one can map the cosmological bouncing problem from the 
Einstein frame to the Jordan frame. The present work shows that a proper understanding 
of the conformal transformations yields interesting properties of the cosmological bounces 
in the Jordan frame. Unlike the case of quadratic gravity, in higher order polynomial 
gravity there are multiple ways in which cosmological bounces can occur. In the case of 
quadratic gravity only bounce in presence of hydrodynamic matter was allowed where 
as in cubic gravity one can also have cosmological bounces in absence of any form of 
hydrodynamic matter. 

The issue of energy conditions was briefly touched in the present work. We have 
shown that one can always produce a cosmological bounce in the Jordan frame without 
violating the weak energy condition and the null energy condition in the Einstein frame. 
The strong energy condition and the dominant energy conditions may be violated near 
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the bouncing point in the Einstein frame. The interesting thing about the Einstein frame 
analysis, of the bouncing problem, is that all the energy conditions are valid at t = 0, 
which corresponds to the actual bouncing time in the Jordan frame. For a cosmological 
bounce the energy conditions in the Jordan frame are violated but one must have to 
remember that the energy conditions are not that clearly and uniquely formulated in the 
Jordan frame. As a result one may not be very serious about the energy conditions in the 
Jordan frame. 

In this work we have introduced the techniques to understand the origin of multi¬ 
ple scalar potentials in the Einstein frame for cubic or higher order f{R) gravity. We 
have explicitly worked out the cosmological bounces arising out of the various potential 
branches in the Einstein frame. Our discussion also includes the approximate Einstein 
frame potential which one can use to study cosmological behavior near i? = 0 as given in 
[32|. As soon one includes the multiple potentials in the Einstein frame the richness of the 
cosmological scenario becomes apparent. In the penultimate section we have presented 
the ansatz following which one may map from the Einstein frame picture to the Jordan 
frame picture without any confusion. 

The other important observation which is reported in the present work is related to 
the cosmological behavior of the universe as it evolves from a bounce to a state where 
the Ricci scalar attains a value R = Rc where Rc is the root of the equation f"{R) = 0. 
Here the primes denote derivatives with respect to the Ricci scalar. It was noticed that 
the conformal transformations between the Jordan frame and the Einstein frame breaks 
down precisely at the point R = Rc for cubic gravity. In the present work we show that 
the state of the universe when R = Rc can never be connected to the bouncing universe 
in a hnite time. In other words the flat FLRW universe, described by cubic f{R) theory 
in the Jordan frame, would take an inhnite time to attain R = Rc after the cosmological 
bounce. The roots of f"{R) = 0 bifurcates the region of cosmological existence. We 
have presented our results for cubic gravity but we hope that our results are general 
in nature and can hold true for higher order polynomial gravity. This particular issue 
requires further research. More over the points R = Rc are interesting, as out of the the 
two branches emanating from this point in Fig. [T| the one on the right of R = Rc has 
/"(R) > 0 and the one on the left of R = Rc has f”{R) < 0. On the left branch one may 
expect Dolgov-Kawasaki like instability in the low curvature limit of the theory. But as 
this branch also includes the point R = 0 one has the freedom for a cross-over to more 
conventional GR dynamics as the Ricci scalar becomes small in magnitude. Ideally the 
Dolgov-Kawasaki instability cannot be applied to the very early universe, in /(R) theory, 
as at those times the curvature may be very high and consequently one can always use 
the left branch to model cosmological bounces happening in the very early universe. 

Lastly, we have utilized all the conceptual understanding about cubic gravity to ac¬ 
tually solve the cosmological bounce problem. The results are given in Fig. Fig. |^and 
Fig. 1^ These results are obtained by solving the problem in the Einstein frame and then 
mapping the solutions back to the Jordan frame. These solutions show various bounce 
possibilities in the Jordan frame for the case of cubic gravity. We expect the possibilities of 
various kind of cosmological bounces will increase with the order of polynomial gravity as 
more and more Einstein frame scalar potentials will emerge in the conformally connected 
Einstein frame. In this paper we do not present the theory of cosmological perturbations 
as our main aim was to analyze the complicated and rich background cosmological evolu¬ 
tion in polynomial /(R) theories. We will present the perturbations on this background 
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in a forthcoming publication. Although, the use of the Einstein frame as an auxiliary 
frame to calculate the background bouncing dynamics was purely a technical choice it 
turned out that this technical choice yields much valued information about the actual 
Jordan frame cosmological evolution if one can disentangle the multi valued correspon¬ 
dence between the conformal frames. The technique evolved in the present work can be 
valuable as one can use pure general relativistic techniques to solve a problem in a higher 
derivative f{R) theory of gravity. Although the methods discussed in this article were 
particularly aimed to tackle the problem of bouncing cosmologies in the Jordan frame, we 
hrmly believe that many of the techniques developed in this article can also be applied to 
understand the general nature of f{R) theories. 
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